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■ The rate at which galaxy halos form is thought to play a key role in explaining many 
CN , observable cosmological phenomena such as the initial epoch at which luminous matter 

forms and the distribution of active galaxies. Here we show how Press-Schechter theory!! 
J> ■ can be used to provide a simple, completely analytic model of the halo formation rate. 

This model shows good agreement with both Monte-Carlo and N-body simulation results. 
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1 Introduction to Press Schechter Theory 
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Following standard Press-Schechter (PS) theory, we convolve a homogeneous, isotropic Gaussian 
Q-i random field representing primordial density fluctuations with a sharp k-space filter. The filter 
size is related to mass of halo by M = 3/ATT 2 pR 3 @. The overdensity of the filtered field at any 
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point traces a Brownian random walk with the variance of the filtered field being the 'time' 
axisE E§. Regions are considered to have collapsed to form dark halos of mass > M if this 
overdensity exceeds a critical overdensity, 5 c (t). The distribution of halo masses is given by 
the distribution of up crossings of the line 5 = S c (t), a solution of the diffusion equation with 
absorbing boundary condition!!!. This leads to a formula for the comoving number density of 
objects of mass M given time t, n(M,t), the standard PS result : 



n{M, t)MdM = pP(M\8 c (t))dM = 2p exp ' ' V< ' ' 
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We normalise the power spectrum by setting the variance of the density field, filtered with a 
top-hat filter of radius 8/i _1 Mpc, a 8 = 0.64. The critical density used with the top-hat filter is 
that predicted from uniform spherical collapse!, that used with a sharp k-space filter is set so 
that we predict the same number density of halos with both niters at the normalisation mass. 



2 The Halo Formation Rate 

We now wish to calculate the distribution of cosmic epochs at which halos of a given mass 
M are formed, P(t\M). Because of the nature of the random walks if we make S c bound, 
$min < &c < 8 max , it can be thought of as a random variable with uniform probability density. 
We can remove these bounds later without affecting the result. Using Bayes' theorem and 
P(a 2 \5 c ), calculated in the previous section, we obtain the conditional probability P(5 c \a 2 M ). 



Taking the limit as 5 m i n — > and 5 max — > oo and changing the variables to mass M and time 
t, we obtain the probability that a halo of mass M formed in the time interval (t,t + dt) : 
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The predictions of equation |2| are compared with the results of both Monte-Carlo analysis and 
an N-body simulation in Fig. |l| Excellent agreement is found between all results. 




Figure 1: Plot comparing Monte-Carlo (open symbols) and N-body (solid symbols) results 
with PS predictions (solid line) of the formation rate of halos of mass 1.3 x 1O 13 M . The 
Monte-Carlo analysis involved 10 4 halo mergers recorded in ~ 10 6 walks, each consisting of 2 11 
uniform steps in a 2 M . The Hydra N-body, hydrodynamics code@ was used to run a simulation 
with 128 3 dark matter particles in a flat A = universe with h = 0.5 and standard CDM power 
spectrum with shape parameter T = 0.25. Groups of particles were found using a standard 
friends of friends (FOF) algorithm with linking length b = 0.2. 



3 Conclusions 

We have shown that we can use Press Schechter theory to produce a simple analytic model of 
the formation rate of dark halos. Our model uses statistical analysis of the trajectories and 
overcomes the problem of relating a trajectory through the hierarchical build up of structure 
to a particular halo. We have shown that the results of this model agree well with those from 
both Monte-Carlo analysis and from N-body simulations. We are grateful for the use of the 
Hydra N-body code || kindly provided by the Hydra consortium. 
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